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We demonstrate experimentally quantum-inspired, spectral-domain intensity optical coherence 
tomography. We show that the technique allows for both axial resolution improvement and disper¬ 
sion cancellation compared to conventional optical coherence tomography. The method does not 
involve scanning and it works with classical light sources and standard photodetectors. The mea¬ 
surements are in excellent agreement with the theoretical predictions. We also propose an approach 
that enables the elimination of potential artifacts arising from multiple interfaces. 
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Optical coherence tomography (OCT) is a powerful, 
three-dimensional (3D) imaging technique which may be 
operated in both the spectral and the temporal domain 
mm, and it is widely employed for biological in vitro 
and in vivo imaging [3]. The axial resolution of conven¬ 
tional OCT is limited by the spectral bandwidth of the 
light source and by the dispersion of the optical compo¬ 
nents and/or the sample under the test. Because of dis¬ 
persion, especially in fiber-based setups, increasing the 
spectral bandwidth of the light source may not necessar¬ 
ily lead to resolution improvement [3] . To circumvent this 
problem several approaches have been proposed, includ¬ 
ing both numerical mm and experimental techniques 
[THS]. In particular, methods such as quantum-optical 
coherence tomography (QOCT) or chirped pulse inter¬ 
ferometry have been shown to possess built-in dispersion 
cancelation and resolution enhancement; however, these 
techniques generally require sophisticated light sources 
such as single-photon sources or chirped ultrashort pulses 
and advanced detection techniques [HHIS]. Further, they 
operate in the time domain requiring in-depth scanning, 
which results in slow measurement times. 

Intensity-based optical coherence tomography, inspired 
by QOCT but using any classical broadband light source, 
was recently put forward in the time domain [HKin]. The 
method produces improved resolution but is hampered by 
the lack of ultrafast detectors capable of recording rapid 
intensity variations characteristic of incoherent sources. 
However, spectral intensity optical coherence tomogra¬ 
phy (SIOCT) provides a much simpler and cost-effective 
alternative that operates, without moving parts, in the 
spectral domain and can use traditional broadband light 
source of any state of temporal coherence and standard 
detectors [min]. Furthermore, SIOCT requires only a 
minor modification of the conventional spectral-domain 
OCT imaging setup while still providing resolution im¬ 
provement and, more importantly, all even-order disper¬ 
sion cancellation. Yet, no experimental demonstration 
of SIOCT has been reported. In this letter, we respond 
directly to this need and show experimentally the func¬ 
tioning of SIOCT with a classical incoherent source, con¬ 
firming both the theoretical resolution enhancement and 
the (group-velocity) dispersion cancellation. The mea¬ 


surements are in excellent agreement with the theoreti¬ 
cal predictions. Our results open up new perspectives for 
high-resolution imaging using classical light sources. 

The difference between the principles of the operation 
of conventional spectral-domain OCT and SIOCT is il¬ 
lustrated in Fig. [l] Light from an incoherent broadband 
source is divided between the two arms of an interfer¬ 
ometer. One arm serves as the reference while the other 
contains the sample to be characterized. In traditional 
OCT, the spectral intensity resulting from the interfer¬ 
ence of the light in the two arms, recorded at the output 
of the beam splitter, is of the form 

-foCTa(w) ^ |As(a;) -|-iAH(w)|^, (1) 

^ocTb(‘^) + ^r(w)P, (2) 

where oj is angular frequency. Eg and Ar are the com¬ 
plex spectral amplitudes of the electric fields emerging 
from the sample and reference arms, respectively, and the 
subscripts a and b denote at which output of the beam 
splitter the spectral OCT measurement is recorded. The 
phase difference between the reference and sample arm 
fields is caused by difference Az in distances that the 
light propagates in air and the propagation through the 



Dispersive sample 


FIG. 1. (color online) OCT and SIOCT concepts. BS - 50:50 
non-polarizing beam splitter. OCT corresponds to a single 
spectral interference measurement at Spectrometer 1 (or 2). 
SIOCT corresponds to a simultaneous measurement of spec¬ 
tral interferences at Spectrometers 1 and 2. 
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dispersive sample of path length L. Both distances mark 
the total path, including wave transits before and after 
reflection. The intensity of the recorded interferogram 
can then be written as 

/ocTa(cc) = + i|^ (3) 

where S{(jj) represents the spectrum of the light source, 
r = |r| exp(i(^s) is the (complex) amplitude reflection co¬ 
efficient of the sample, and c the speed of light in vacuum. 
Further, j3{u}) denotes the frequency-dependent propaga¬ 
tion constant within the sample. Using a Taylor-series 
expansion for the propagation constant around the cen¬ 
tral frequency loq of the source spectrum up to second 
order, one obtains 

-^OCTa(w) = S{u1q + Uj') 

X -I-i|^, (4) 

where lu = ujq + lu' and r = Az/c + /3iT, with /3o, /3i, 
and P 2 being the propagation constant, the group de¬ 
lay, and the group-velocity dispersion at frequency Wg, 
respectively. The envelope of the spectral interference 
pattern represents the source spectrum whilst the fre¬ 
quency, phase, and amplitude of the modulation under¬ 
neath depend on the sample position and reflectance, 
which can then be obtained through a Fourier transform 
of the recorded interferogram. The resolution is given by 
the full width at half maximum (FWHM) of the Fourier 
transform of the term corresponding to the optical path 
difference between the reference and sample arms. In 
the absence of dispersion, the axial resolution is inversely 
proportional to the source bandwidth. With dispersion 
present in the system, the resolution decreases by a factor 
[1 -I- (/32T)^(Aa;/2-\/2In2)^]^/^, where Aw is the FWHM 
spectral bandwidth of the light source. 

In SIOCT, the spectral interference patterns of the 
fields in the two arms are recorded simultaneously by 
two separate detectors at the two output ports of the 
beam splitter, and the location and the reflectance of 
the sample are obtained from the Fourier transform of 
the cross product of the individual spectral intensities 
C{lo') = /ocTa(wo +w')-foCTb(wo - w'). It is straightfor¬ 
ward to show that the interference pattern of the SIOCT 
signal is then given by 

= S'(a;o -I- uj')S(uJo — w')[co(a;') -I- ci(w') -I- 02 ( 0 ;')], 

( 5 ) 

where 

co(w')= (|rp + l)=^ 

-l- } (6) 

ci(w') = -4(|r|2 + l)3{e-‘™'} 

X }, (7) 

and 

C2(a;') = -2|r|25i{e-'2™'}, (8) 


with 5ft and 3 denoting the real and the imaginary parts, 
respectively. One sees that the SIOCT signal produces 
an interferogram with an envelope bandwidth reduced 
by a factor of \/2 compared to that of the conventional 
OCT inteferogram envelope. The Fourier spectrum of the 
function C(lij') gives access to the sample information and 
consists of three separated peaks corresponding to the 
Fourier transforms of co(w'), ci(a;'), and 02(0;') convolved 
by the Fourier transform of S (wg -I- w' ) S' (wg — w' ). Several 
observations can be made: 

We first remark that the product S(a;g-|-a;')S(wg—w') is 
always an even function, even if the spectrum of the light 
source is not symmetrical with respect to wg. The term 
cg(w') is an even function whose Fourier transform pro¬ 
duces a real-valued peak centered at the zero delay (equal 
path lengths). The term ci(a;'), on the other hand, is an 
odd function that depends on the optical path length and 
its Fourier transform corresponds to an imaginary-valued 
peak centered at the optical path difference between the 
two arms. Finally, the Fourier transform of C 2 (w') gives 
rise to a real, negative-valued peak at twice the optical 
path difference. It is precisely this peak in the Fourier 
spectrum that gives information about the sample posi¬ 
tion and, after correcting for the distance, one obtains 
an overall resolution improvement of as compared to 
standard OCT. This also means that the imaging depth 
is only half of the one obtained by the standard OCT, 
as the number of data points is constant and their den¬ 
sity is doubled in SIOCT. Note further that in practice 
the terms arising from Ci(a;') and 02 ( 0 ;') can be distin¬ 
guished as one is real and the other imaginary. Besides 
the resolution enhancement, another benefit of SIOCT 
is the inherent dispersion cancellation of all even-order 
terms in the Taylor-series expansion. This is because in 
SIOCT intensities of opposite frequencies relative to cen¬ 
tral frequency wg are multiplied, canceling the even-order 
phase terms that arise from dispersion. 


Light source Imaging spectrograph 



FIG. 2. (color online) Modified Hong-Ou-Mandel interfer¬ 
ometer, experimental setup: BS - 50:50 non-polarizing beam 
splitter. SI, S2 - spectrograph input and output slits. Dl, 
D2 - photodiodes. 

We confirmed experimentally the resolution improve- 
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merit and dispersion cancellation as compared to stan¬ 
dard OCT using a modified Hong-Ou-Mandel (HOM) 
interferometer m, as shown in Fig. With our setup, 
the signal recorded at either detector independently is 
identical to a conventional OCT system, while measur¬ 
ing simultaneously the signal at both detectors allows to 
construct the SIOCT interferogram. The light source is a 
fiber-coupled superluminescent diode (Exalos ESL 1620) 
with a center wavelength of about 1610 nm. The source 
spectrum is close to Gaussian with a spectral bandwidth 
of 55 nm (FWHM). The light emitted by the source was 
collimated with a parabolic mirror and three identical 
non-polarizing beam-splitter cubes were used to divide 
the light among the sample and reference arms. A sin¬ 
gle, partially reflecting mirror placed on a manual trans¬ 
lation stage to adjust the optical path length was used 
as the sample. The reference wave transmitted through 
beam splitter BSl was reflected from an additional beam 
splitter BS2 in order to equalize the dispersion and phase 
shifts experienced by both beams due to the various op¬ 
tical elements. As a result, any difference between the 
complex amplitudes of the electric helds in the sample 
and reference arms before interfering on beam splitter 
BS3 is caused only by the path difference between the 
two arms and sample presence. Light at the two output 
ports of BS3 was coupled to a single-mode hber using 
achromatic lenses. The hber outputs were then placed on 
top of each other in the object plane of a monochroma¬ 
tor (Horiba iHR 550) allowing to measure simultaneously 
the spectral intensities at the two detectors. On passing 
through the output slit of the monochromator, light was 
collected by multimode hbers and intensities were mea¬ 
sured by InGaAs amplihed diodes (Thorlabs PDAIOCS). 
Lock-in detection was used to improve the signal-to-noise 
ratio. Wavelength scanning and data acquisition were 
controlled by a PC. 

In order to conhrm the resolution improvement we hrst 
performed measurements in the absence of dispersion in 
the sample arm. The results are shown in Fig. where 
both the real and the imaginary part of the Fourier trans¬ 
form of the cross-product function S{uj') are presented. 
For comparison, the theoretical results obtained on the 
basis of Eqs. (|4) and ([^ are superimposed as circles. The 
optical path de ay was converted into physical distance as 
light would travel in vacuum and thereby the 0 point rep¬ 
resents equal path lengths. For ease of comparison with 
the conventional OCT result, the distance was divided 
by 2 for the SIOCT measurement so that the image peak 
would correspond to the actual sample position. In gen¬ 
eral we observe excellent agreement between the exper¬ 
imental and theoretically predicted results. Specifically, 
we see that the measured imaginary part of the Fourier 
spectrum is an odd function and corresponds precisely to 
the Fourier transform of the term ci(a;'), while the mea¬ 
sured real part is even and matches closely the Fourier 
transforms of the terms co(a;') and 02 ( 0 ;') of the SIOCT 
interferogram. Significantly, we also observe how the res¬ 
olution is improved in comparison to the standard OCT 


measurement (see inset in Fig. [^. The FWHMs of the 
OCT and SIOCT peaks corresponding to the sample po¬ 
sition, summarized in Table |lj are 23.5 /rm and 17.3 /rm, 
respectively, which gives a ratio of 0.736, close to the 
theoretical resolution improvement predicted for a light 
source with a Gaussian spectrum. 



FIG. 3. (color online) Comparison of conventional OCT and 
intensity-based SIOCT in the absence of dispersion in the 
system. The solid lines show the experimental results and the 
circles correspond to theoretical predictions, (a) OCT result, 
(b) and (c) Real and imaginary part of the Fourier transform 
of the SIOCT interferogram, respectively. The inset in (b) 
displays a direct comparison of the axial resolution of both 
techniques: black OCT, blue SIOCT. The intensity units have 
been normalized in this case. 

We next proceeded to confirm experimentally the in¬ 
herent dispersion cancellation of the SIOCT scheme. For 
this purpose, an 8 cm thick bulk piece of SFIO glass was 
inserted into the sample arm. In order to increase the 
effect of dispersion the beam was made to propagate 8 
times through the glass cube, effectively corresponding 
to a 64 cm glass piece. The group-velocity dispersion 
coefficient of SFIO at 1610 nm is f }2 = 1.63410“^®s^/m, 
which should lead to a resolution decrease by a factor 
of about 2.4 in the case of conventional OCT. The mea¬ 
surement results for conventional OCT and SIOCT are 
illustrated in Fig. As for the dispersionless case, the 
theoretical results computed from Eqs. @ and are 
also superimposed as circles. It can readily be seen that 
the experimental results in both cases follow the theoreti¬ 
cal predictions with great accuracy. More specifically, we 
observe how in the case of conventional OCT the FWHM 
of the Fourier transform peak corresponding to the sam¬ 
ple location has broadened nearly 2.5 times to 58.0 /rm 
due to the dispersion of the SFIO glass, a value very close 
to that expected from the theoretical calculations. Most 
remarkably, the width of the SIOCT Fourier transform 
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FIG. 4. (color online) Comparison of conventional OCT and 
intensity-based SIOCT when dispersion is present in the sys¬ 
tem. The solid lines show the experimental results and the cir¬ 
cles give the theoretical predictions that are based on Eqs. 0 
and 1^. (a) OCT result, (b) and (c) Real and imaginary part 
of the Fourier transform of the SIOCT interferogram, respec¬ 
tively. The inset in (b) illustrates the direct comparison of 
the axial resolution in the two techniques: black OCT, blue 
SIOCT. The intensity units have been normalized in this case. 


peak representing the sample position is unaffected by 
the presence of the bulk piece of glass, clearly demon¬ 
strating that dispersion is inherently canceled in SIOCT. 
In fact, the resolution was found to be 16.2 ^m, which 
is even slightly better than in the absence of dispersion. 
This slight change in the resolution is caused by third- 
order dispersion, which is not canceled in SIOCT. The 
measured FWHMs, 16.2 /rm and 58.0 /im, of respectively 
the OCT and SIOCT peaks corresponding to the sample 
position are summarized in Table |lj together with the 
dispersion-free case. Finally, we also see how the artefact 
peak corresponding to the term ci(w') is generally af¬ 
fected by dispersion, which provides an additional means 
of discrimination. 

In the case of multiple interfaces, however, the arti¬ 
facts have generally complex values due to cross-talk be¬ 
tween the terms corresponding to the multiple interfer¬ 
ences. This means that they can not be discriminated 
from the real interfaces by solely considering the real 


part of the Fourier transform of the SIOCT interfero¬ 
gram. However, in this case, one can take advantage of 
the fact that the amplitude of the artifacts correspond¬ 
ing to the ci(a;') term in Eq. ^ depends on the center 
frequency wq of the source spectrum. Indeed, for a sta¬ 
tionary light source the spectrum need not to be sym¬ 
metrical with respect to ojq, so that one can choose arbi¬ 
trarily the center frequency as long as the detected signal 
significantly exceeds the measurement noise. Because the 
phase of the ci(uj') term oscillates with wq, potential arti¬ 
facts will then vanish when multiple measurements per¬ 
formed with spectra of different center frequencies are 
averaged. Interestingly, multiple measurements with a 
different center frequency can be numerically performed 
by post-selection of a smaller number of points from the 
measured interferogram, which in turn shifts the spectral 
window together with the central frequency. In this way, 
one can produce an ensemble of data sets corresponding 
to spectra with different central frequencies and averag¬ 
ing the Fourier spectra eliminates any artifact due to the 
presence of multiple interfaces. Yet, it is important to 
emphasize that averaging has its limitations as the dis¬ 
persion coefficients Pi depend on the frequency ojq around 
which the Taylor expansion is performed. This leads to 
the situation where different central frequencies give dif¬ 
ferent sample positions and, in that case, averaging may 
lead to broadening of the image peak and a decrease in 
the overall resolution. 


TABLE I. Summary of resolution comparison 


SIOCT 

Standard OCT 

Without dispersion 17.3 ^m 
With dispersion 16.2 ^m 

23.5 m 

58.0 /rm 


In summary, we have experimentally demonstrated 
spectral intensity optical coherence tomography and 
proved the associated resolution improvement and dis¬ 
persion cancellation. We also suggested a possible ap¬ 
proach to eliminate artifacts that might be present in the 
case of samples with multiple distinct interfaces. No fast 
detectors are needed in our technique as only mean spec¬ 
tral intensities are measured. The method works with 
classical light sources and standard photodiodes, illus¬ 
trating the simplicity of the technique and showing its 
potential application to high resolution imaging. 
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